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Abstract 

The collective behavior of molecular motor proteins have been investigated in the litera- 
ture using models to describe the long-time dynamics of a unidimensional continuum motor 
distribution. Here, we consider the phenomena related to the transport of particles (vesicles, 
organelles, virus, etc) in the realm of these continuum motor systems. We argue that cargo 
movement may result from its ability to perturb the existing motor distribution and to surf at 
the resulting shock waves separating regions of different motor densities within the transient 
regime. In this case, the observed bidirectionality of cargo movement is naturally associ- 
ated with reversals of shocks directions. Comparison of the quantitative results predicted 
by this model with available data for cargo velocity allows us to suggest that geometrical 
characteristics of the transported particle shall determine the extension and intensity of the 
perturbation it produces and thus, its dynamics. Possible implications of these ideas to virus 
movement at the cell body are discussed in connection with their distinguished morphological 
characteristics. 

Key words: collective effects of molecular motors; cargo transport; non- linear partial dif- 
ferential equations; shock waves. 



1 Introduction 



The intracellular active transport of particles, including organelles, vesicles and virus is mediated 
by motor proteins such as myosin, kinesins and dyneins [1] . The unidirectional motion of a single 
motor along protein filaments is well characterized experimentally and it was first modeled at a 
microscopic level by a stochastic dynamics describing the behavior of a Brownian particle in the 
presence of a time-dependent asymmetric potential field [2] , [3] , [4] , [5] . The idea is based on the 
mechanism of "ratchet and paw " introduced by Feynmann to discuss the meaning of the second 
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law of thermodynamics [6]. Since then, this model has been used as a prototype to explain why 
and in what conditions Brownian particles are able to do work against external potential gradients. 
Also if particles can bind directly to motors, the transport of cargo follows as a direct consequence 
of motor movement. 

More recently, it appeared in the literature attempts to describe the movement of interacting 
Brownian motors, since it was realized that collective effects emerging from this situation may 
be relevant to explain certain characteristics of cellular transport and, in particular, the observed 
property of bidirectionallity of cargo movement [7], [8], [9], [10]. This non-diffusive type of process 
is characterized by inversions of cargo direction after processive runs 1 that may be preceded by 
relative large resting times (intervals within which the particles remain at the same microtubule 
site) . 

Attempts to describe the observed bidirectional movement of extensive objects (filaments) me- 
diated by motor proteins use essentially the same ideas of the original microscopic models [11]. 
Numerical results indicate that some of the observed characteristics of bidirectional movement are 
captured in this way provided that the number of motors attached simultaneously to the filament, 
arranged in a periodic way, be no less than a certain critical number N m i n . Because of this, it is 
not clear how such models could be extended to describe the observed bidirectional movement of 
small particles, as vesicles or virus. 

The most accepted explanation in these cases is referred in the literature as the coordination 
model [7] . According to this, the bidirectional movement would result from the coordinated 
action of two type of motors - a plus-ended and a minus-ended motor - attached simultaneously 
to cargo (filament). The reversal of its direction would just reflect the fact that one or other type 
of motor, but not both simultaneously, shall be active during the respective time intervals. The 
motor coordination, that is, the control of motors activity would be accomplished by an external 
non-motor protein complex that should be able to coordinate the movement and timing of many 
motors of different nature and distinct characteristics - certainly a very non-trivial job. Such an 
external complex has never been identified in any in vivo experiment [12]. 

At certain scales of interest, phenomena related to the collective behavior of interacting Brow- 
nian motors can also be approached by models that are intended to describe dynamic aspects of 
continuum media. Such description, referred in the following as macroscopic model, is generally 
justified upon evaluation of the characteristic sizes and time scales where the molecules operate at 
very low Reynold's numbers [13]. Usually, the approach is based on the continuum versions of the 
"asymmetric exclusion processes" (ASEP) [14], [15], [16], [17] for studying the long-time behavior of 
motors interacting through volume excluded [13] [18] [19] [20]. From this perspective, the individual 
microscopic asymmetric movement of motors is assumed a priori and for open boundary conditions 
the problem consists in analyzing the steady-state behavior of a defined motor density as the so- 
lutions to the corresponding non-linear differential equation - the non-viscous Burgers equation - 
that describes the dynamics of the system in these limits. It is also possible to superimpose to the 
ASEP a non-conservative Langmuir process to allow the system exchange motors with the bulk at 
any position of the microtubule [19]. Apparently, however, questions related to cargo transport, 
have not been considered in this context . 

Here, we make a proposal in this direction based on considerations about the interactions 

1 Processive run in this context is referred to the movement accomplished in one definite direction before inversion 
or detachment from filament. 
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between cargo and motor molecules. It relies on the idea that the arrival of a cargo particle may 
perturb the motor system in such a way to produce local changes onto the existing motor density. 

Thereafter cargo may take advantage of the gradients induced by this initial perturbation to 
move along microtubule by surfing at the density shock waves formed as the motor system relaxes 
back to the situation before interaction with cargo. Shock waves separate regions of different motor 
densities and evolve according to the transient solutions to the corresponding Burgers equation for 
the considered initial conditions (perturbation). Within this view, bidirectional movement of cargo 
particles follows naturally as a consequence of the reversals of shocks propagation direction. 

From a more quantitative perspective, we use data from the movement of vesicles in Drosophila 
embryo [21] to perform a phenomenological analysis of the model considering the explicit expressions 
obtained for shock velocities. As we shall see, this allows us to associate geometrical characteristics 
of cargo with the extension and magnitude of the perturbation it produces. We then discuss on 
the possibility to make a connection between these results with the properties of the observed 
movement of virus at the cell body. 

2 From the microscopic potential models to the macro- 
scopic motor density profile 

We consider the continuum limit of the mean field approximation (the macroscopic limit) of an 
ASEP discrete model representing the stochastic dynamics of particle motors moving in a one- 
dimensional space. Within this view, the diffusion process is not accounted for explicitly; its effects 
on the movement of each motor under the time- dependent asymmetric potential, are assumed a 
priori. Yet, the parameters that define the asymmetry of the potential, can be incorporated into 
the parameters of the considered ASEP. 

Consider the time-dependent asymmetric and periodic microscopic potential shown in Fig. 1. 
The position of each minimum may be associated to the position of a one-dimensional lattice site; 
the distance between sites being I = a + b coincides with the spacial period of the potential. The 
difference between parameters a and b define the strength of the asymmetry. The time dependence 
is such that it alternates between an "off" state (flat potential), when the particle feels no forces, 
and an state "on" - the sawtooth profile shown. If the potential is turned off and the particle 
starts to diffuse from a local minimum corresponding to the site i of the lattice, it shall diffuse 
around this position up to the instant when the potential is turned on again, after a time interval 
t. At this instant, the particle is supposed to fall immediately into the nearest minimum (adiabatic 
approximation). Because a > b, and considering a regime of relative slow diffusion, it is more likely 
that the nearest minimum achieved by the particle be that at site i + This prototype mechanism 
was conceived to explain from a microscopic point of view the observed unidirectional movement 
of a Brownian motor [3], [4], [5]. 

ASEP models have been introduced more recently in this context as a way to account for collec- 
tive effects of many motors moving on the same lattice model ([18], [20]). At these scales, motors 
are self-driven and the interactions among them are that of volume excluded, i.e. a motor can not 
occupy a site already occupied by another motor. Let us now consider a stochastic description of 
the process at this scale. The time unit being the interval r within which a particle can execute 
one movement to a neighbor site with (transition) probability pi to move from site i to site i + 1 



3 



and (transition) probability to move from site % to site i — 1. Neglecting correlations (mean- field 
approximation) the average density piit + r) of motors at site i, at time t + r satisfies the following 
recurrence relation: 

Pi (t + r) = Pi (t) + p i+1 (t)q i+1 (t) + - Pi(t)(pi(t) + qi(t)) (1) 

Excluded volume interaction is introduced into the model by attributing to pi and qi an explicit 
dependence on the occupation of the target sites: 

Pi (t)=p(l-p i+1 (t)) (2) 

and 

qi {t) = q(l-p i . 1 (t)) (3) 

Here, p and q are dimensionless constants such that p, q < 1. Now, let N be the total number of 
lattice sites and let L = 1 be its total length. As in Ref.[19] , the discrete model can be coarse- 
graining with lattice constant I = 1/N to a continuum such that for iV — * oo the density pi(t) 
becomes a function p(x,t) of a continuum variable x = i/N. In this limit, pi±i can be expanded in 
powers of I : 

p i±1 (t) = p[x ±l,t)= p{x, t) ± W x p(x, t) + h 2 d 2 x p{x, t) + 0(Z 3 ) (4) 
and the recurrence in (1) becomes 



p(x,t + r) -p(x,t) 

= 7 



T 

where we have defined 



W 2 x p(x,t) + (p-q)(2p(x,t) - l)d x p(x,t) 



+ 0(/ 2 ) (5) 



7 = - (6) 



For r — > and Z — > keeping 7 in (6) finite and neglecting terms of 0(/), Eq. (5) converges to 
the non-viscous Burgers equation 

d t p + 1 (p-q)(l-2p)d x p = 0. (7) 

Now, a connection between the above description and the microscopic model can be made by 
setting 

p — a/l and q — b/l (8) 

so that Eq. (7) is rewritten as 

d t p + Kd x (p(l-p)) = Q (9) 

where we have defined K = 7(0 — 6)iV positive. This can also be expressed as d t p + d x j = with 
the particle current given by 

j(x,t) = Kp(x,t)(l-p(x,t)). (10) 

Notice that to keep K finite as iV — > 00, parameters a and 6 must scale as 1/N. The same behavior 
is expected for the Langmuir rates in the model of Ref. [19]. In these conditions, and with no lost 
of generality, we can set K — 1. 
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The nonequilibrium steady-state solutions to this equation have already been explored in the 
references mention above to study the long - time behavior of molecular motor traffic with open 
boundary conditions. This allows one to make important predictions on the stationary density 
profiles related to the concentrations gradients of motors at the microtubule ends. In particular, 
the emergence of stable domain walls separating regions of different motor densities, predicted in 
these conditions, have indeed been observed in recent experiments [22]. Here, we consider periodic 
boundary and study the short-time behavior of the system described by Eq. (9), after being 
perturbed with respect to the motor density steady-state p(x,t) = 1/2. The results we obtain for 
a particular choice of initial conditions are shown in the next section. 



3 Results 

It is known that perturbations around the steady state solutions of Burgers equation (9) may create 
shocks due to encounters of particles at different velocities ([23], [24]). This leads us to speculate 
that cargo may use these shock waves, by surfing at the wave fronts, in order to reach target sites. 

To illustrate these ideas, we can generally conceive that the approach of cargo induces a pertur- 
bation on motor density such as to create at the cargo neighborhood regions where motor density 
becomes lower and others where density becomes higher than the steady state p = 1/2. As an 
example, we choose 



p(x,0) 


= 1/2 


for 


x < 


p(x,0) = 


1/2 -e 


for 


< x < 2a 


p(x,0) 


= 1/2 


for 


2a < x < 3a 


p(x, 0) = 


1/2 + e 


for 


3a < x < 4a 


p(x,0) 


= 1/2 


for 


x > 4a 



to represent such perturbation (see Fig(2a)). Parameters a and e are related, respectively, to the 
extension and magnitude of this perturbation. 

We then use the method of characteristics to find the transient solutions to Eq. (9) that reveal 
the time evolution of shocks for these initial conditions (IC) (11), which is our main interest here. 
This is performed in the Appendix. 

The results for the density profile are represented in Fig. (2) at various instants of time. The 
behavior of each shock front, namely, distance traveled up to an eventual encounter with other 
shock fronts - which are identified here with the processive runs - the time spent within each of 
these runs, and the corresponding expressions for average shock velocities are compiled in Tablel. 

In the next Section, we use these results to examine the consequences of the hypothesis made 
above. 



4 Discussion 



5 



4.1 the surfing 



The model proposed here offers rather simple explanations to a few but essential phenomena related 
to motor-driven cargo transport. In particular, it suggests an alternative to interpret data showing 
that cargo executes a non-diffusive bidirectional movement before it arrives at its final destination. 

In general the main difficulty to explain the bidirectional movement of cargo relies precisely on 
the fact that the existing models have been conceived to explain the unidirectional movement of 
motors. Because of this, the original idea that the cargo movement shall be put into effect by its 
attachment to a single or a few identical motor proteins has been reviewed in the literature. 

The alternative mechanism proposed here does not exclude the possibility that eventually, 
the presence of two types of motors may contribute to explain the reversals of particle direction. 
However, it does not use this as a necessary condition as does the coordination model mentioned 
above. Within such an alternative context, cargo is driven by shocks and this is independent of the 
nature of motors involved provided they offer an assessable density background. That is, in order 
to move and eventually change direction, cargo must be able to perturb this background and follow 
gradients of motor density right after perturbation, just like the putative microscopic "second class 
particles" introduced in the context of many-particle asymmetric exclusion processes (ASEP) to 
study the nature of shocks at these scales [14], [15]. 

Although we do not offer here any direct evidence that would support these ideas we can test 
for their compatibility by confronting some experimental data with the results shown in Table 1 for 
the considered initial conditions. For example, it is shown there that a cargo particle accompanying 
the wave front referred here as wave front 2, or second shock, would perform a backward movement 
within a time interval Ti = a/2e that depends on the model parameters. Within the subsequent 
time interval T 2 = a/4e, the particle would remain at rest up to the instant 3a/ 4e when it reverses 
the direction of movement. Thereafter, the particle would proceed eventually changing the mag- 
nitude of its velocity though not direction anymore (in this specific case (11) the reversal happens 
just once) until the motor density recovers its steady state profile p(x,t) = 1/2 and the particle 
reaches its final destination at the one of the microtubule ends, as expected. Therefore, at least 
qualitatively, the phenomenological aspects of the observed movement, as reversals of direction, 
processivity in both directions and resting times are accounted for by this model. 

Evidently, if one considers other initial conditions, for example by introducing additional pa- 
rameters to characterize the perturbation, it is expected a great variety of time intervals and 
corresponding distances traveled in both directions, differing in processivity and velocity. There- 
fore, a relevant aspect of this proposal is the possibility to attribute primarily to cargo the task 
for creating such gradients, thus inducing its own surfing movement. This is discussed below in 
connection to a more quantitative analysis of available data from real biological systems. 

4.2 phenomenological analysis 

Consider for instance the data in the literature on the movement of lipid vesicles on Drosophila 
embryo [21]. We can check for their compatibility with the results in Table 1 all expressed in terms 
of just two parameters a and e. To accomplish this, we observe that data analysis obtained from the 
transport of vesicles in the referred experiment identifies two typical runs in both directions: one 
that is short (average distance d short ~ lOOnm) and slow (average velocity v s i ow ~ 200 nm/s) and 
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another that is long (average distance d\ ong ~ lOOOnm) and fast (average velocity Vf ast ~ 400 nm/s) 
. We shall use these data to estimate the values of a and e. We proceed by identifying v s i ow , taken 
from experimental data, with the average velocity of the 4 th shock front v± = 2 e/7 . It results 2 
e ~ 700. We can now use this to predict the magnitude of the average shock front velocity for the 
5 th shock t> 5 = e/2 ~ 350 which approaches very well the values (absolute values) observed for Vf ast 
~ 400nm/s. 

On the other hand, these results must be compatible with the predictions for the distances trav- 
eled by the corresponding shocks. This complements the analysis by estimating possible values for 
parameter a. Accordingly, the predicted distance must be approximately one order of magnitude 
greater than d±. This is in fact accomplished by the results shown in Table 1 since from there 
d^/d^ = 8. 

We shall now perform a slight different analysis of the data by redefining the constant a such 
that 

a = 500p(nm) (12) 

where p is a numerical factor. From this, d± = 80p(nm) and d$ = 650p(nm). Therefore, if we choose 
p = 1.2, it results d^ ~ 100(nm) and d^ ~ 800(nm). These can be compared to the results for d s h or t 
and di on g given above to conclude that the characteristic times and distances traveled by shock 
fronts 4 and 5 can be related in a very direct way with the typical values obtained in experiments 
for short and long runs, respectively. 

More interesting, however, is the fact that this approach to the experimental data is made 
possible with the choice in (12) for p ~ 0(1). First, notice that 500nm is the average diameter of the 
vesicles transported in the referred experiments. Thus, if the nature of interactions between cargo 
and motors are of short range, this implies that the magnitude and extension of the perturbation 
should be related to the geometrical characteristics (typical sizes) of the cargo. Therefore, if the 
nature of interaction of different cargos with motors is the same, but cargo differ from each other 
in morphological aspects, these differences would be reflected in the specific way (size, magnitude) 
each of them perturb the equilibrium motor distribution and consequently, in the diversity of 
movements. Such diversity is indeed observed in experiments using different motor/cargo systems 
and/ or under different conditions [7]. 

The above considerations seem specially attractive to explain the behavior of virus particles 
moving within the cellular environment observed more recently using confocal microscopy tech- 
nique ([25], [26], [27]). According to the authors of these experiments, virus appear surfing along 
filopodia in such a way to perform a bidirectional motor-driven process in the presence of a single 
type of motor protein, unidirectional myosin II. From these observations the authors then discuss 
on the possibility that, in these experiments, virus might have been able to create the conditions to 
induce their own movement. These observations support the idea that virus use (hijack) existing 
cell transport machinery to invade the cell [27]. The model presented here suggests a mechanism 
to explain how this can be accomplished without external and specific mechanisms of coordination. 
Such explanation has a focus on the morphological aspects of cargo - which, of course, is a distin- 
guished characteristic of virus. In our view, this complements the suggestions made in Ref. [25] 

2 It must be remembered that the parameter e was defined adimensional and that the velocities are expressed in 
units of the parameter K that have been defined equal to one. 
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relative to the existence of a common mechanism that drives different virus particles to their targets 
in cells, in spite of their morphological differences. We propose that this common mechanism is 
just related to their ability to interact with the system by perturbing the existing motor profile. 
Specific geometries, however, would be used by virus to create a variety of movements, possible 
necessary for their efficiency We thus believe that these ideas may be helpful for studying the 
mechanisms of virus infection and also for designing drugs that must be directed to the cell body. 



5 Appendix 

Here we obtain the properties of shock waves shown in Tablel within the transient regime, after 
the motor system is perturbed by interaction with cargo. We apply the method of characteristics 
to find the transient solutions to the equation (9). For this, we follow the texts by Habermann 
([23]) and that by Evans ([24]). 

+ (1-2^,0)^ = (13) 

for the chosen initial conditions (IC) (11). In turn, these solutions define the properties of shocks 
depicted in Table 1 which is our main concern in this work. For each pair (p(xo),xo) there corre- 
sponds a characteristic curve for (9) given by 

x(t) = (l-2p(x ))t + x (14) 

Thus, the IC in (11) define the different characteristic families: 

x < (15) 



Xi(t) 


= x 


for 


x 2 (t) 


= Xq + let 


x 3 (t) 


= Xq 


for 


x 4 (t) 


= x 


- 2et 


x 5 (t) 


= Xq 


for 



4a < Xq 

These are represented in Fig. (3) where one observes the occurrence of shocks dividing regions of 
different density values. There are six shocks. There are also shown in this figure the rarefaction 
regions that are localized between any two characteristic families that get apart from each other as 
the time evolves. The exact determination of the origin of each of these shocks (time and space) - 
and rarefaction regions and the course of their time evolution allows us to describe the evolution 
of p(x, t). 

At the initial time (t = 0) there are two shocks: 

• shock 1 (x s i(t)) at position x sl (0) = 2a, between the families defined by x 2 (t)and x 3 (t). The 
velocity of this shock front is determined by mass conservation conditions : 

dx sl j(p(x 3 )) -j(p(x 2 )) 



dt p(x 3 ) - p(x 2 ) 



e 
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where the current j(p) is defined by (10). We then have 

x sl {t) = 2a + et (16) 
that gives the shock position as a function of time. 

• shock 2 (x s2 (t)) at x s2 (0) = 3a, between the families defined by x 3 (t) and x A (t). The velocity 
of this shock is given by 

dx s2 = j(p(x 4 )) -j(p(x 3 )) = _^ 
dt p(x A ) - p{x 3 ) 

from which 

x S 2(t) = 3a — et (17) 
Still at t — 0, there originate two rarefaction regions 

• rarefaction 1 (x r i(t)) at x r i(0) = 0, evolving in the region between families xi(t) and 
x 2 (t). The requirement of entropic solutions drive us to look for the corresponding family of 
characteristic curves x r i(t) defined by the density p r \ at each point in this region: 

x rl (t) = (1 - 2p rl )t for 1/2 -e < p rl < 1/2 (18) 

which are also represented in Fig. 2. Conversely, one can write the solutions to p r i{x r i,t) for 
density profile at each position x r \ in this region as 

p rl (x rl ,t) = l/2-— for 0<x rl <2et (19) 

• rarefaction 2 (x r 2(t)) at x r 2(0) = 4a, evolving in the region between families Xi{t) and 
xs(t). The characteristics in this region are given by 

x r2 (t) = 4a + (1 - 2p r2 )t for 1/2 < p r2 < 1/2 + e (20) 

thus 

p r2 {t,x r2 ) = 1/2 + — — — for 4a - 2et < x r2 < 4a (21) 

• shock 3 (x s3 (t)) The third shock occurs at time t 3 when the two families x sl (t) and x s2 (t) 
cross. Thus, t° is determined by the condition:a; s i(t3) = ^2(^3), which gives 

t° 3 = a/2e. (22) 

The velocity of the shock is given by 

dx S 3 = j(pM) ~ j(p(x 2 )) 
dt p(x 4 ) - p(x 2 ) 

meaning that the shock front is stationary at x s3 (t 3 ) = ^1(^3) = 5a/2. 



= 
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• shock 4 (x s4: (t)) occurs at time t° A when the family x 4 with origin at x = 4a reaches the 
second region of rarefaction and also cross the characteristics from family x<i- Thus, t\ is 
determined by the relation x s3 (t1) = x r2 (p = 1/2 + e; t®), which gives 

t\ = 3a/4e (23) 
The (continuity) condition, in this case reads — = j_^L_!£_^ — i'(p( 4)) ^ ^ 

dt Pr2{XsAj ~ P{Xa) 

dx si 4a-x s4 

~df = e — — (24) 

then x S 4 is given by the solutions to the above differential equation for initial condition x S 4,(t1) 
such that x s s(t1) = x S 4(t^). The result is [28] 

x s4 (t) = 4a - 2v / 3aet + 2et (25) 
The dependence on time of the corresponding shock velocity can now be found 



^ = 28- J 3 -^ (26) 
dt V t 

• shock 5 (x s5 (t)), occurs between the two rarefaction regions, at the time t®, when the 
characteristic from the family X2 that originates at xo = encounters x S 4 

x r2 (pri = 1/2 - e; tg) = x s4 (t° 5 ) (27) 

resulting 

n 4a , , 

4 = S (28) 

The mass conservation conditions in this case ( see 21 and 19) reads 

dx sb _ j(pr2(Xs5)) -j(Prl(x s5 )) _ X s5 - 2d 
dt Pr2(x sb ) - p r l(x s5 ) t 

Considering that x S 4(tf) = x s5 (t®) this can be integrated to find 

x s5 (t) = 2a + {e/2)t (29) 

thus 

IF = 5/2 (30) 

• shock 6 (x S 6(t)), occurs at ig when shock 5 encounters the characteristics from family x$ 
that starts at x = 4a. Then, t$ results from the relational (ig) = 4a. Using the above derived 
expression for x s5 (t) one gets 

t° 6 = Aa/e (31) 
The mass conservation conditions impose that 

dXs6 = jjpM) -j(Prl(Xs5)) = ^£6 

dt p(x 5 ) - p rl (x s5 ) 2t 
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whose solution for the initial condition x s6 (tg) = 4a becomes 

x s6 (t) = 2Vaei (32) 
from which one can find the shock velocity: 

dx S Q las 



dt 



(33) 



The fact that this velocity is always positive, implies that for t > = Aa/e the solutions to 
Eq. (9) for the considered initial conditions (11) are 

p(x,t) = 1/2 for x<0 (34) 

x / as 

p(x,t) = 1/2-— for 0<x<4a + J— (35) 

2£ y t 

/ CLE 

p{x, t) = 1/2 for 4a + J — < x (36) 

From these, one sees that for long times, i.e. for t > t$ = 4a/ e, the density profile returns to its 
uniform steady state profile, p(x,t) = 1/2. This corresponds to the situation before perturbation 
imposed by cargo particles. The solutions discussed above refer to behavior of p(x, t) within the 
transient regime soon after being perturbed by cargo. The properties of shocks derived here are 
compiled in Table 1. 
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Table 1 
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Table 1 - Kinematic properties of shock fronts formed within the course of time evolution of the 
initial density profile in (11). The average velocity is calculated from the ratio between traveled 
distance (column 4) and the corresponding time interval (column 3) for each processive movement. 
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6 Figure Caption 



Figure 1 - The asymmetrical potential field introduced in Refs. ( [3], [4], [5]). 
This is defined by parameters a and b and we chose a > b. 

Figure 2 - The motor density profile at different instants of time. The cargo is represented by 
the circle and at each instant, its position is adjacent to the shock front, the arrows pointing to the 
direction of the movement. Fig.2(a) shows that the considered initial condition (11) corresponds 
to a local perturbation on motor density with respect to its steady state p(x,t) = 1/2. 

Figure 3 - Characteristics to Eq.(9) for the initial conditions given by (11). 
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